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OO ' Abstract 

f^i , An effective action is obtained for the area and mass aspect of a thin shell of radiating self- 

0^ ' gravitating matter. On following a mini-superspace approach, the geometry of the embedding 

space-time is not dynamical but fixed to be either Minkowski or Schwarzschild inside the 
shell and Vaidya in the external space filled with radiation. The Euler-Lagrange equations of 
motion are discussed and shown to entail the expected invariance of the effective Lagrangian 



^ | under time-reparametrization. They are equivalent to the usual junction equations and 

I , suggest a macroscopic quasi-static thermodynamic description. 

5h . 



1 Introduction 



The Einstein field equations for a general distribution of matter are a formidable challenge and 
the study of the collapse of self-gravitating compact objects is a very hard task to which much 
effort has been dedicated since the birth of General Relativity (see, e.g., Ref. Q). One of the 
first papers in this field showed that, when gravity overcomes all other forces, a sphere of 
matter collapses under its own weight into a point-like singularity. This opened up a whole line 
of investigation about the nature of such a singularity and the way it forms. 

A state corresponding to a point-like singularity would violate Heisenberg's uncertainty prin- 
ciple, therefore such a problem makes physical sense in a region for which gravity can be treated 
at the (semi) classical level (see e.g. [||, f|, || and Refs. therein). Provided the conditions for 
such an approximation hold, one is left with two types of difficulties which conspire against 
the achievement of a definitive answer: firstly, one would like to consider a realistic model for 
the collapsing matter and, secondly, one has to face the intrinsic non-linearity of the Einstein 
field equations. The former aspect includes the quantum nature of matter as described by the 
standard model of elementary particles with all of its own intricacies. The latter difficulty be- 
comes greater along with the complications of the matter model but, in turn, provides the most 
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interesting features. Thus one has to find a sensible compromise between reality and practical 
solvability of the equations. 

A great deal of simplification usually follows from global space-time symmetry. The natural 
framework for the study of gravitational collapse is isotropy with respect to a point. Of course 
this rules out rotating objects, such as stars and other astrophysical systems, and freezes many 
of the degrees of freedom of gravity (Birkhoff's theorem forbids the emission of gravitational 
waves Q). However this does not render the field equations trivial and some of the strong field 
effects certainly survive in this approximation (see, e.g., Ref. ||). 

An example of manageable isotropic distributions is given by dust fluid models, such as the 
Tolman-Bondi space-time j?]] , which can be used to represent spheres of pressureless matter and 
cosmological models. Such space-times develop caustics, where the energy density of the fluid 
diverges, because the geodesies followed by dust particles cross for generic initial conditions. 
The dynamics can be investigated by approximating the continuous distribution of matter by a 
discrete medium of nested homogeneous time-like shells of small but finite thickness. When dust 
particles are confined by their own weight to roughly within a few Compton wavelengths and the 
latter is negligible with respect to any other "macroscopic" length in the model (to wit, the square 
root of the area of the shell and its Schwarzschild radius ||] ) , one can take the limit of infinitely 
small thickness {thin shell limit) and treat each shell as a singular hypersurface generated by 
a ^-function term in the matter stress-energy tensor [10|] . Then the basic problems are the 
study the dynamics of one shell and the collisions between two nearby shells (an issue we do not 
touch upon here). 

The method for treating singular hypersurfaces in General Relativity was formulated in 
Ref. H as following from the Einstein field equations. It amounts to the Lanczos junction equa- 
tions between the embedding metrics and entails the conservation of the total stress-energy. The 
equations of motion of thin time-like shells have then been derived for all spherically symmetric 



embedding four-geometries [11], including the vacuum Q and bubbles of non- vanishing cosmo- 
logical constant [12]. It is remarkable and physically sound that, in non-radiating cases, the 
trajectory of the area of the shell is determined by the junction equations once one has fixed 
the two external metrics and the equation of state for the density and the surface tension (the 
radial pressure must vanish). The special case of radiating shells was studied in Ref. |l*|] for 
unpolarized radiation of relatively high frequency which behaves as null dust and gives rise to 
a Vaidya external geometry ||T4|, [Dm. This is indeed a good approximation when the radiation 



wavelength is much smaller than any other "macroscopic" length in the model (see, e.g. [16 
and Refs. therein) and is reasonably consistent within our approach, since the same condition 
on the Compton wavelength of the matter in the shell is required by the thin shell limit ||. 

For the purpose of developing a quantum theory, knowledge of the equations of motion 
is not sufficient, instead one needs an (effective) action for the shell degrees of freedom. The 
conceptual starting point is thus the Einstein-Hilbert action rather than Einstein field equations. 
In Ref. |l(| both metric and matter degrees of freedom were kept dynamical and the general 
form of the action was given for a barotropic fluid with step- or 5-function discontinuity. One 
then implements the symmetry of the system in order to reduce the action to a manageable 
effective form. The literature on this topic treats essentially two approaches for the shell: 

Canonical approach. In Ref. [17], the embedding empty space-time is foliated into spatial 
sections of constant time according to the ADM prescription [18], with lapse and shift 



functions in the four-metric. The Einstein-Hilbert Lagrangian from Ref. [fil]] is then in- 
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tegrated over the spatial sections by making use of the properties of a spherical vacuum 



IE]. This leads to a canonical effective action for the canonical variables of the system 
which include the inner and outer Schwarzschild times on the shell as reminders of the ge- 
ometrodynamics of the embedding space-time. A set of constraints ensures the invariance 
of the action under reparametrization. The shell matter Lagrangian used in jl7]] is not in 
the form of a field theory and does not give rise to any canonical variable. 



Mini- super space approach. In Refs. [|l^, 20 1 the embedding metrics are chosen a priori to 
be specific solutions of Einstein field equations and are not dynamical. The corresponding 
contribution to the Einstein-Hilbert Lagrangian can then be integrated over a convenient 



space-time volume and expressed in terms of metric variables of the shell world- volume [20|. 
This yields an effective action which is invariant under reparametrization of the time on 
the shell and is equivalent to the canonical one in domains of the phase space where some 



of the constraints are solved classically [21]. Thus one expects that the mini-superspace 



approach gives a limited information on the quantum theory of geometry fully addressed 
in the canonical approach, but does not limit the possibility of taking into account the 
quantum nature of matter in a semi-classical context, as was later done in Ref. ||. 

In the present paper we illustrate the derivation of an effective action for a time-like shell 
which can emit unpolarized high frequency radiation. Since the basic aim is to establish a 
starting point for the study of semi-classical effects induced by the quantum evolution of the 
matter in the shell (along the lines of Ref. J|]), the mini-superspace approach will suffice. The 
dynamical meaning of such an action will just be to generate the time evolution of the shell 
degrees of freedom, which will be identified with the area and mass aspect of the shell. However, 
we shall also see that this approach, unlike the standard junction prescriptions, allows a simple 
interpretation of the evolution equations and time-reparametrization invariance in terms of the 
thermodynamics of the shell. 

The plan of the paper is as follows: in Section we start from the general Einstein-Hilbert 



action as given in Ref. [1C] and simplify it in order to derive a mini-superspace effective action. 
In Section || we write the Euler-Lagrange equations of motion, two of which are equivalent 
to the junction equations and compare with previous approaches. In section |] we illustrate 
a thermodynamic interpretation and in Section |5] we finally comment our results. We shall 
follow the sign convention of Ref. JIJ (see also Appendix |b|) with Greek indices /i, u, . . . labeling 
space-time four-coordinates and Latin indices i,j,... for the coordinates of a three-dimensional 
sub-manifold; k = 16ttG, with G the Newton constant. 



2 Derivation of the action 

The space-time f2 we are considering is parted into two regions, Q.^, separated by the shell 
world- volume S. inside the shell is devoid of matter and f2 + possibly contains out-flowing 
radiation. The corresponding Einstein-Hilbert action, from Ref. |T(|, is thus 



Seh = 




+- / afxVh[K]t- I d 3 xVhp 

+- I d A x^lZ+ f d 4 x^£ rad , (2.1) 
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Figure 1: The space-time of the collapsing shell with the ADM foliation. 



where g is the determinant of the four-dimensional metric, h the determinant of the three- 
dimensional metric on E, 7^- is the scalar curvature, K the trace of the extrinsic curvature of 
E, p the energy density of the shell and C ra d the radiation Lagrangian density. We have also 
introduced the notation [F]_ for the difference between the limiting (on the shell) values of a 
function F computed in Vt + and f2~. 

The above action is the starting point for our derivation. First of all we observe that the 
spherical symmetry of the system makes it possible to introduce three global coordinates (r, 6, (ft), 
with r > the radial coordinate of a sphere of area Airr 2 , 9 G (0,7r) and 4> G (0, 2ir) the usual 
angular coordinates. The next step is to decide which one of the two approaches outlined in the 
introduction is better suited for the present problem. Since a formulation of the kind given in 



Ref. [19 1 is not available at present for a radiation-filled space-time p^ j, the canonical approach 
is not viable. Therefore, we follow Ref. |20| and assume the Einstein field equations are satisfied 
inside the space-time volumes [23]. 



The above assumption fixes the form of the embedding metrics. In particular, we observe 
that Q can be naturally parted into three regions (see Fig. |l|): 

i) the inner empty space Qj n = whose geometry can be either Minkowski or Schwarzschild. 
Let ti n be the (Schwarzschild) time coordinate in f2j n , then the corresponding metric can be 
written 

ds 2 n = -a in dtl + a" 1 dr 2 + r 2 dtf , (2.2) 

where cftl 2 = d6 2 + sin 2 9 d(j) 2 is the usual line element on the two-sphere and aj n = 1 — 2 Mj n /r. 
The constant Mj n is positive for Schwarzschild and zero for Minkowski; 

ii) an outer space £l ou t which can be filled with radiation. If t ou t denotes the (Schwarzschild) 
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time in Q ou t, the metric, which is Vaidya |l4j], takes the diagonal form 

ds 2 mt = b 2 out a~l t dt 2 out + a^ t dr 2 + r 2 dQ 2 , (2.3) 

where b ou t = dtm/d r m, a ou t = 1 — 2m/r. The mass aspect, m = m(t out ,r), equals the total 
energy included inside the sphere of area 4 tt r 2 at time t out Q ; 

Hi) an empty space fioo prior to any emission of radiation with Schwarzschild metric and mass 
parameter equal to the total ADM mass of the system, 

ds 2 ^ = -aoo dt 2 ^ + dr 2 + r 2 dtt 2 . (2.4) 

where Oqo = 1 — 2 M^jr. 

In this frame, the shell trajectory is given by r = R(t) and t in / out = T in / out (r), where r is 
the arbitrary time variable on the shell world-volume S with three-metric h{j given by 

dsl = -N 2 dT 2 + R 2 dn 2 , (2.5) 

N = N(t) being the shell lapse function. The full stress-energy tensor of the matter in the 
system contains two parts, 

T„ v = 5(r -R) + T;f 9{r - R) , (2.6) 

where 5 is the Dirac 5- function and 9 the step-function. The source term, 

Sij = diag \n 2 p, -R 2 P, -R 2 sin 2 6 p\ , (2.7) 

is the three-dimensional stress-energy tensor of a fluid with density p and (surface) tension P 
and T™ d is the stress-energy tensor of the out-flowing null radiation. 

The mass aspect can be extended to all values of r > 0, does not decrease for increasing r 
and must be continuous in f2 except at the shell radius, thus 

lim m(T in ,r) = M in 

lim m(T out , r) = M out > M in 

r— >R+ 

lim m(too, r) = M oc > M out , (2.8) 

! >00 

where M out = M out (r) is the mass aspect at the (outer) shell radius and equals the total energy 
of the shell plus Mj n . Other matching conditions will be implemented in the following and, as 
in Eq. (|2.8j ), we shall use capital letters for the restriction (or limit) to the shell position of 
functions of the space-time coordinates which are denoted by the corresponding small letters, 
e.g., Ai n = ai n (R,T). We shall usually drop the subscripts in/out whenever this does not cause 
any confusion (e.g., t = ti n for r < R) and total derivatives with respect to r are denoted by a 
dot. 

Since the metrics above solve the Einstein equations with the source (|2.6| ), the volume con- 
tributions in the action (2.1) must vanish identically. Further, the matching between £l ou t and 
Sloo is smooth, provided the radial coordinate R s of the border between the two regions satisfies 



dR s 2 Moo , s 

dT = 1 " ' (2 ' 9) 
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where equals M ou t at the time r at which the emission starts, and no boundary terms arise 
at the surface r = R s . The reduced Einstein-Hilbert action can then be written 

S%* = — HdrNR 2 [/CC*-47r H drNR 2 p + S B 

= S s G heU + Sff 11 + S B , (2.10) 

where Sq 1 ^ 1 is the surface gravitational action of the shell, Sb contains all surface contributions 
at the border of the space-time £1 (including those which are usually introduced to cancel second 
derivatives of the dynamical variables) and Sf^ el1 is the shell matter action, 

gshell = _ 4?r [ Tf drNR 2 P , (2.11) 



where we assume p does not depend on N and its dependence on the other two shell degrees of 
freedom R and M out will be clarified in Section |3[ 

At this point one can foresee a technical difficulty in evaluating surface terms at t constant 
because the region fl ou t is not homogeneous (see Fig. ||). In fact, the mass aspect m is constant 
along the lines of the outgoing flow of radiation and one can introduce a null coordinate u = 
Uout{t ut-,f) such that out-going null geodesies are given by u constant with four- velocity ~ 
(0, 1, 0, 0) and one has 

T rad --- — — S u <) u (2 12) 

^ " nr 2 du » v ■ { ' 

The mass aspect m = m(u) does not depend on r now, and the metric is written (for r > R) 

ds 2 out = -a out du 2 -2dudr + r 2 dti 2 . (2.13) 

Such a u ou t is defined also for M ou t constant, in which case it coincides with a null Eddington- 
Finkelstein coordinate for the Schwarzschild space-time [|]J (see also Appendix |A|). Integration 
over r at time t ou t constant would then involve integrating functions of m and its derivatives, 
that is the knowledge of M(T) and R(T), for all t{ <T < t ou t- However, although the differential 
relation between (t out ,r) and (u ou t,r) is known, 

du = -a~ut (bout dt out + dr) , (2.14) 

u = u ou t(t ou t,r) cannot be written explicitly unless the functions M out (r) and R(t) are fixed 
once and for all. At the same time the trajectory of the shell is to be derived from the sought 
effective action and cannot be given a priori. This renders explicit integration over slices of 
constant t inconvenient. 

The above argument suggests one define a "mixed" foliation, where volume sections are 
defined by t constant for r < R and u constant for r > R (see Appendix [A] for an equivalent 
overall null foliation). Hence we shall integrate the Einstein-Hilbert action in the four- volume 
= VLi n U flout displayed in Fig. Such volume has a time-like boundary at r = i?^ and 
two "mixed" boundaries defined by U n = U, u ou t = Ui and ti n = tf > t{, u ou t = Uf- The shell 
trajectory is described equivalently by r = R(t), t = T(r), with (fixed) end-points Ri = i?(rj) 
at T(ji) = ti and Rf = R(tj) at T(tj) = tf, in £7j n and r = R(t), u = U(t) in flout, with 
end-points at U(ti) = Ui and U(jf) = Uf. 
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Figure 2: The space-time U £l ou t with the "mixed" foliation used for the computation of the 
effective action. 

The fact that all volume terms (the integrals over Jl 1 * 1 in Eq. (^j])) do not contribute to the 



reduced action (|2.10[) can now be checked explicitly. The basic observation is to note that the 
scalar curvature for a metric of the form ( 2.1 3| ) or ( |A.1| ) is given by (see Appendix [B]) 



d ( r 2 da 
dr \ 2 dr 



and vanishes identically once one substitutes in a 
The corresponding volume action is thus 



47T 
K 



tf 



dt 



R 2 , „ 4 ^ 

r dr TZ in H 



a ou t for Vaidya or a 



2 

du I r drlZout 
IR 



(2.15) 



for Schwarzschild. 



, 



(2.16) 



as expected. One also recalls that the Lagrangian density for a field of null dust is proportional 
to kfj, k^, where is the fluid four- velocity, and thus vanishes along null geodesies. 

In order to compute Sq we introduce Gaussian normal coordinates (r, 77, 9, </>) near the 
shell such that the trajectory R(t) is given by 77 = with 77 > for r > R. Then, the jump in 
the components of the extrinsic curvature of the shell world-volume are given by the relation |l| 



dgij 



dr/ 



1 %j 

2 dr/ 



(2.17) 
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In detail (see Appendix |C|) 



sin 



p Pin 

~~N 



R 



P 



(2.18) 



out 



r\ = +e 



JC 7 



where 



Pi; 



■ . iv 3 a^ in 



RN-^^-RN 
2 <9i2 



2 I AwiV 



&4 



out 



(2.19) 



r] = +e 



The above expressions give 
8vr 



nshell 



(It 



2(5R + 



P 



p= \JAN 2 + R 2 



■ N N 2 
R - R N + - 



(2.20) 



dA p-R \ dA 
dR + \ AN j dU 



(2.211 



out 



where we used dai n /dui n = to write the limits from the two regions in the same form. 

The integral in Eq. ( ^2l|) contains a second derivative of R which should be eliminated by a 
suitable term in Sb- Such term must be a total derivative so that it does not affect the equations 



of motion and for the Schwarzschild metric it can be found in Ref. [12]. The same kind of term 
works for Vaidya as well and is given by 



7<1) 



8 7T f T f d d 
K J T , dr 



R tanh 



-i 



Adding fl2~2p to ( ^2l|) gives 



-ishell 
'G 



+ s 



(1) 



8vr 



d,T 



2PR-2RR tanh 



-i 



r; 
p, 



R 

P. 



(2.22) 



out 



+ 



MP MR 



A 



A 



(2.23) 



- out 



We observe that Pi n /Ai n = Ti n [|T^] and the mass aspect Mj n is constant, therefore the last two 
terms are irrelevant in J7j n (being zero or total derivatives). In O, ou t, however, they cannot be 
neglected. 

We shall now identify another surface term of dynamical relevance which arises from the 
borders of O. We observe that the expression ( p. 15 ) of the scalar curvature is defined up to the 
addition of zero which can be written in the form 



2a d . . 2jdfr 2 da 

= — — — (r — ar) + 



dr 



- 2 dr \ 2 dr 



(2.24) 
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where a and 7 are arbitrary constants. Hence 

r 2 dr 1Z = 



4ir 


f u f 


rRoo 




/ du 




K 




Jr 


8?r 


f u f 


rRoo 




/ du 




K 




Jr 



dr 



;i + a)-(l + a) Q-(ar) 



d ( r 2 da 



(2.25) 



from which one readily sees that there are no (non-trivial) surface terms at the null borders 
u = Ui and u = Uf of OU £. It is shown in Ref. fl2j that no such terms arise along the spatial 
borders t = tj and t = tf of fij n . Therefore, from (|2.25[) one concludes that a surface term could 
come only from the border at r = i?oo and should be compensated by 



,(2) 



8vr 
K 

K 



du 



(1 + a) ar + (1 - 7) 



r 2 da 
~2 dr~ 



du [(1 + 2 a + 7) M OU i - (1 + a) i? c 



(2.26) 



For a = 1 and 7 = one obtains the usual contribution (see Appendix |D[). However we find it 
more useful to set a = — 1 and 7 = 0, for which the integrand in Eq. ( |2.26| ) with the relation 
( C.6 ) becomes 



S 



(2) 
B 



8vr 

K 



r f Pout 

dr 



R 



A, 



M, 



O il I 



out 



(2.27) 



since this form cancels one of the terms in the integral ( |2.23| ). 

Putting together all the relevant pieces given in Eqs. ( |2.11| ), ( |2.23| ) and ( 2.2 7| ) finally yields 
the effective action 



S eff (N,R,M t 



out J 



2Rf3 -IRRt&nhT 1 f- 



+ 



M mi , R M nut R 



- out 



A 



out 



A 



out 



-4vr [ Tf dTNR 2 p 



drL eff {N,R,M out ) , 



(2.28) 



for the variables N, R and M OU f It is worth noting that one recovers the effective action for 
a shell which collapses without emission of radiation simply by setting M out = 0. In fact the 
integral of M ou t R/A ou t, which does not appear in the action given in Ref. f20|| , depends only on 
the end-points when M out is a constant and is then dynamically irrelevant. 

In the next Section we shall vary S e ff in order to check that it generates the same equations 
of motion as following from the junction conditions reviewed in appendix |E|. 



3 Equations of motion 

In the following we shall obtain the Euler-Lagrange equations of motion by varying the effective 
action with respect to the r-dependent quantities appearing in the effective Lagrangian with 
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the end-points held fixed at T{ and t/. Such variations are formally computed by acting on the 
Lagrangian with the operator (i = 1,2,3) 



d d ( d 



8X i dXi dr \dX* 



(3.1) 



with X i = (N, R, M out ) and S&fc) = SX^Tf) = 0. 

To begin with, we note that it is sensible in the spirit of General Relativity to assume that 
the shell proper energy E does not depend on the shell velocity R. Further, in order to keep a 
sufficient degree of generality, we consider E depending on both shell variables R and M out (but 
not on N) , so that the physical shell energy is related to the quantity 4 tt R 2 p by 



E = 4irR 2 p{R,M ( 



out ) , 



(3.2) 

and agrees with the source in the 00-} unction equation ( |E.2[) . 

We next introduce the canonical momenta Pj = (dL e jf/dX l ), which are formally given by 



Pn 
Pr 



Pm — — 




2R 



tanh 



R 



+ 



M, 



out 



out 



A 



out 



K A, 



(3.3) 
(3.4) 

(3.5) 



out 



The particular expressions ( |3.3|) and ( |3.5D deeply affect the nature of the dynamical system at 
hand, as can be seen from the (symmetric) matrix 



Wi. 



dL 



dP 



dX { dXi dXi 



(3.6) 



whose rank equals the number of canonical degrees of freedom with Dim[VF] — Rank[PF] being 
the number of primary constraints. Since 



RankfVF] = Rank 



diag , 



dP R 
OR 







1 



(3.7) 



the Lagrangian L e ff contains two primary constraints and is said to be non-standard [24 |. 
Varying S e ff with respect to the lapse function yields 



5L 



eff 



dL 



eff 



16 TT R 



5N 



dN 



K N 



]m 
out 



E 



(3. 



where we made use of our previous assumption that p does not depend on N. Upon setting this 
variation to zero we obtain the first primary (Hamiltonian) constraint 



H 

N 



E 



16 TT R 



K N 



m 
out 



o , 



(3.9) 



which is related to the time-reparametrization invariance of the shell three-metric and is formally 
equal to the analogous constraint obtained in Ref. (2^] for non-radiating shells. In order to make 
contact with the notation in Appendix |E] we observe that, for N = 1, one formally recovers the 
09-} unction equation ( E.2 ). 
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Varying with respect to R yields the equation of motion 

SL ef f _ dL eff d (dL e ff 



5R 



dR 



dr V dR 



. 



that is 
5 

Jr 



(NE) 



16 7T 



P + 



ft 



RR + 



N 2 M RNR 



R 



N 



+ 



M, 



out 



(3, 



out 



R 



out 



Pout l-2M out /R 



Upon defining 



P = 



1 



JL(MF)- — 
\itNR~5R^ >~Ja 



(3.10) 



(3.11) 



(3.12) 



(A = 4:irR 2 is the area of the surface of the shell) and setting TV = 1 we formally recover the 
rr-j unction equation (|E.3|) . 

Now, unlike the non-radiating case, one must also consider varying M ou t, thus obtaining the 
second primary constraint which has no interpretation as a junction condition, 

SL eff _ dLeff d / dLeff 

dM dr \ dM 
16 7r (3 out - R 5 



SM, 



out 



l-2M out /R 5M, 



out 



{NE) = 



(3.13) 



It was already obvious from Eq. ( |3.3[) and W\\ = that the lapse function N is a Lagrange 
multiplier and can be assigned any function of r. Hence, from now on we work in the proper 
time gauge N = 1, which makes most of the equations look simpler. The corresponding con- 
straint ([O]) was in anticipation named the Hamiltonian constraint since the quantity N H is the 
canonical Hamiltonian, as can be seen from the Hamiltonian form of the effective Lagrangian, 

L eff = P N N + P R R + P M M out -NH. (3.14) 

The Hamiltonian constraint must be preserved in time, which yields the secondary constraint 



dH 
~dT 



dE_ 

AttR 2 



dE 
~dR 



R 



dE 



dM, 



M, 



out 



out 



2 — (p — ) — — ^ out P° ut ~ 

— { - p) - - R R _ 2 M ^ 



, 



(3.15) 



where we have used both Eqs. ( |3.11 ) and (3.13) to show that, in general, Eq. ( |3.15 ) is trivially 
satisfied by our ansatz E = E(R, M out ) in Eq. (|3.2|) . 

Since W33 = 0, also the quantity M out is not a true dynamical variable. Further, the 
total time derivative of ( 3.13| ) vanishes in virtue of E = E(R, M ou t) on using Eqs. ( 3.13j ) itself 
and (3.15), therefore no new secondary constraint is generated. This means that one can set 
M ou t = Moutir) (any function of r), which determines E = E(R, M ou t) through the constraints 
(|3.9| ) and (|3.13j) , and compute R = R(t) with correspondingly Rfa) = Ri by solving Eq. ( 3.11 
This picture is well suited to describe a shell whose luminosity curve in time is given by [|l3| 



dQ 
dr 



16 vr ft 



'out 



R dM, 



out 



l-2M out /R dr 



(3.16) 
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The above quantity is negative for radiating shells (we required U > 0, see Eq. ( C.6 )) and 
diverges for the shell approaching its Schwarzschild radius ]25|], unless 

dM out 



M, 



out 



u 



du 



{R-2M ( 



out I 



(3.17) 



for R — > 2M out . This singular behaviour (and the corresponding bound on M out ) is due to 
the fact that the coordinates (u, r) are related to the point of view of a static observer (see 
Appendix |e|) and is removed by passing to Israel's coordinates [|l5|] (see also Appendix |C|). 

However, one can alternatively choose the value of E = E(R(t), M out {r)) as an explicit 
function of r and determine P from ( 3.11| ) and the luminosity from ( |3.13j ). Then the trajectory 
R = R(t) is obtained by imposing the Hamiltonian constraint ( |3.9| ) for all times Tj < r < tj. 
This option is particularly useful if one wishes to impose only initial conditions and then consider 
both the trajectory and the luminosity of the shell as completely determined by the specific 
interaction between the shell matter and the emitted radiation encoded in the dependence of E 
on R and M out . 



3.1 Comparison with other approaches 

In a previous approach [jl2] , ^0f the outer mass is considered fixed, in which case it is not possible 
nor necessary to consider the variation with respect to it. Nontheless one may proceed in analogy 
with such an approach and obtain a self-consistent, complete set of equations starting from our 



action Eq. (|l2Sj) . 

The first case one may consider is that for which the arbitrary function M ou t(r) is taken to 
be a function of the shell radius M ou t = M ou t(R(r)). One then has that the two terms to the 



right of eq. ( 2.28 ) are a total derivative with respect to the proper time, 



(It 



M out R M out R 



A 



out 



A 



out 



(- 



log A 



out 



(3.18) 



which is a boundary term and does not influence the equation of motion. Then the tension is 
found to be 



p= 2 - 

K 



PR 1 dA 
R + ]3 + 2~pdR 



(3.19) 



- out 



which is the analogue of eq. ( |3.11| ) and agrees with the TT-junction equation. 

Using this expression and imposing the conservation of the Hamiltonian constraint we obtain 
as anticipated a continuity equation, which takes the form 



dR R y 



P) 



(3.20) 



Finally for the general case M ou t(r) we note that Eq. ( 3.13| ), which was obtained by varying 
the action with respect to M ou t, can also be arrived at by using the expression for the pres- 
sure (that is the equation of motion for R) and imposing the conservation of the Hamiltonian 
constraint 

dH 

- = 0. (3.21) 



It is clear, as we have mentioned after Eq. ( 2.28 ), that for M out constant we reproduce the usual 
results HH]]. 
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4 Thermodynamics 



Just as for the black hole case [26], one may attempt for the case of the radiating shell a 
thermodynamic approach. The first observation is that in thermodynamics one has "quasi- 
static" processes so that the equations of state are always satisfied. This of course implies in 
our case that all time derivatives be considered small so that matter on the shell always be in 
equilibrium J27|. Let us illustrate such an approach (in this Section we denote the quasi-static 
limit of previously defined quantities by the same letter in italics style). 

In the shell one expects that the intensive coordinate be the surface tension V and the 
extensive one be the area A. Thus, considering quasi-static processes for which 



~o 2 IVlnitf 

« 1 ■ 



from Eq. ( 3.11| ) one finds (again setting N = 1) 



V 



2 



R 



M 



R\fA. 



out 



Similarly, from Eq. (E 



one may identify 
„ 16tt 



R 



out 



(4.22) 



(4.23) 



(4.24) 



which will correspond to the shell internal energy. We note that all quantities are a function 
of the (slowly) time-dependent variables M ov ± and R, thus it is clear that d£ (a change of the 
internal energy) is an exact differential [28] of M out and R, which is a statement of the first law 
of thermodynamics 



One may now introduce a variation of the heat Q, 



d£ 
dS 



TdA 

8rcRVdR 



— — dM out + — 
dM out \dR 



8irRV) dR 



which on using Eqs. (14.23|) and flj~2l) leads to 

16 7T dM out 



dQ 



y/l-2M out /R 



(4.25) 



(4.26) 



This is in agreement with the quasi-static limit of Eq. (|3.13j ) and is related to the luminosity of 
the shell, Eq. ( [3.16 ). At this point, in order to obtain the change in entropy, one must introduce 
a temperature T which must be such that dS = & Q/T is an exact differential as a statement of 
the second law of thermodynamics [p9|. This implies that 



which has a solution 



_d_ 

dR 



T 



TJl-2M out /R 



%Kk b M out y/1-2 M out /R 



(4.27) 



(4.28) 
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where kj, is the Boltzmann constant, and corresponds to the desired equation of state connecting 
Mout, R and T. The entropy change in this case is then given by 

dS = (^-^-) 8 vr M ut dM out , (4.29) 



V ti k 

and for k = 16 tt, h = = 1 one has (omitting an arbitrary constant) 

S = 4 vr M^ ut = - (area of horizon) , (4.30) 

which is the usual result for a black hole and it is clear that a collapse for which M ou t = is 
adiabatic ( "isentropic" ) . 

Obviously, given our expressions for the entropy, surface tension and internal energy, one 
may actually evaluate the various thermodynamic potentials. Let us just end by determining 
the specific heat at constant area. It will be given by 

which for R > 3 M out is negative as expected. 



5 Conclusions 

In this paper we have analyzed the dynamics of radiating shells in General Relativity by deriving 
a mini-superspace effective action for the area and mass aspect of the shell. Besides proving 
the dynamical equivalence of this approach with the usual treatment via junction equations, we 
have introduced a temperature (equation of state) and obtained a thermodynamic, quasi-static 
description for the evolution of the shell. 

The effective action in itself is useful for developing the quantum theory and we hope to use 
it along the lines of Ref. (8|. We also think that the remarks in the Introduction and the analysis 
of the equations of motion carried on in the previous Section make clear that description given 
for radiating shells is a general framework which can be used to study a wide variety of cases of 
physical interest. 



A Overall null foliation 

Instead of the "mixed" foliation used in the text, one might introduce a null Eddington- 
Finkelstein coordinate [l| diii n = dti n — a" 1 dr such that 

dsj n = —ai n du in — 2 dui n dr + r 2 d£l 2 , (A.l) 

for r < R, and consider slices of constant u both inside and outside the shell. Then the volume 
over which one integrates the Einstein-Hilbert action would have null boundaries at u = U{ and 
u = Uf for all < r < R^ (this would change the shape of f2j n with respect to Fig. It is 
however easy to show that this does not change the effective action. In fact the scalar curvature 
for Schwarzschild and Minkowski can be formally written as in Eq. ( 2,15| ), with aj n a function 



of (at most) r, since the mass aspect is constant and equal to Mj n for r < R. Therefore both 
the volume contribution and surface terms vanish in the region r < R and one is left with the 



terms at r] = — e in as given in Eq. ( [2.2 
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B Curvature scalar for the Vaidya metric 



The scalar curvature for a metric of the form (|2.13| ) or ( |A.1| ) with a generic a = a(u, r) can be 
computed straightforwardly from the definition of the Riemann tensor (TJ 



where ]i,v . . . = u,r,9,(f), and 



-t->m — pM _ p/ 1 ip" rf _ pA 1 p« 



The non- vanishing connection coefficients, 



1 



are given by 



pu 



1 da 



uu 2 dr 



T u\ = n 9 ^ ^crv,X + 9aX >" ~ 9uX ^ 



IT. = sin 2 6 r&> = r 2 sin 2 i 



-pr 1_ Oa. _j_ a d<i -pr 

1 uu — 2 du 2 dr ' 1 4 



sm 2 en. 



-ar sin 



->r 1 9a 



(B.l) 
(B.2) 

(B.3) 



(B.4) 



p0 _ p^ _ 1 p0 __ : 2«p<? 

1 r6» — 1 r</> — r ' <t><t> ~ b 

and lead to the result stated in Eq. (2,15), 



8111 V COS I 



C Gaussian normal coordinates in Vaidya space-time 

The metric in a neighborhood of the shell can be written in Gaussian normal coordinates 
(t, 77, 0,0) as 



OOO 00 

ds = g TT dr + g m dr] + 2 g rr? dr dr] + r dQ , 



(C.l) 



where 



0tt(V, 77) = -a - 2 



<9u \ 2 ^du dr 
dr J dr dr 

du du du dr du dr 
dr] dr dr dr] dr] dr 
du\ 2 _ du dr 



ar V (T,v) = = 

1 fdu\ z ^du dr 

w ' \drj J dr] dr] ' 



(C.2) 



and we choose 77 > (77 < 0) for points in fi ouf (fij n ). 

Matching ( |C.1| ) to the three-metric fl2,5| ) on the shell world- volume at 77 = gives the set of 
equations 



AU 2 + 2UR = N 2 



dr] 



A 



du 
dr] 



u ^ r 

\ 2 



■ du 



r)=0 



2 ^ U 
77=0 / 



dr 

r,=0 &n 



-1 . 



(C.3) 
(C.4) 

(C.5) 



11=0 
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Prom ( |Cl| ) it follows that U = (±0 - R)/A, with given in Eq. (pHo) ). We require U > and, 
since we wish to describe collapsing trajectories with R < for i? > 2M, we must choose the 
plus sign 



U 



0-R 



(C.6) 



For R — > 2 M the above expression diverges (unless = R = 0), which signals the fact that the 
coordinates in use become singular on that surface p5|]. 

In fact, one has that u — ► +oo for r->2m (both for Vaidya and Schwarzschild). This can 
be cured from the onset by passing to Israel's coordinates (v, w) which are regular all the way 
down to r = [lq|. They are defined by 



1 



W(u) ' dv 4m(v) 
and r = 2m(v) + W(u) to. In this frame the Vaidya line element becomes 



ds 2 



4 dm 



+ 



w 



W dv 2mr 



dv 2 + 2dvdw + r 2 dVL 1 



(C.7) 



(C.8) 



However the same result for the components of the extrinsic curvature is obtained if the change 
is made at the end of the computation (see Appendix ||). 
Eqs. 



and ( C.5 ) yield 

du 
dr] 



r)=0 



0-R 
' AN 



Or 
drj 



T)=0 



N 



(C.9) 



Hence the radial coordinate r as a function of r and rj is continuous, but its derivative with 
respect to r\ has a jump on E. At the same time u need not even be continuous across S. We 
can now write (u, r) as functions of the Gaussian normal coordinates (r, 77) explicitly up to order 



u 



U 



0-R 
AN 





(CIO) 



R + j^V + Otf) 



where we recall that the above expressions hold both in Qi n (rj < 0) and O, ou t {v > 0), thus it is 
to be understood that 



.4 



lim a in = 
lim a out 



1 



2M m /R 
-2M out /R 



in 9,i n 
in £l ut 



{CM) 



and so forth. 

From the knowledge of the mapping ( |C.10 ) one can determine the components ( |C.l| ) of the 
metric and their derivatives with respect to rj and compute the extrinsic curvature according to 
Eq. 
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D Standard boundary term at Roo 

Surface contributions are usually computed by making use of the trace /C of the extrinsic cur- 
vature of the border of the space-time volume according to Eq. fl2.10| ). Further, K, is related 
to the covariant derivative of the unit normal to the border. Thus, at r = Roo, one has 



S 



(3) 
B 



tij 



dudOdd) VhX7„n>* 



R 



(D.l) 



where h = ar A sin 2 9 is the determinant of the pull-back of the metric fl2.13j ) on an hypersurface 
of constant r, V denotes the covariant derivative in the metric ( 2.13j ), 



V u n» 



-9 



-gW 



(D.2) 



and = (—a l l 2 , a 1 / 2 , 0, 0) is the unit normal to an hypersurface of constant r in (u,r) com- 
ponents. Computing the derivative at r = i?oo one finds 



s (3) = 8vr 
B k 



du 



2R n 



dM. 



out 



Rn 



2M, 



out 



du 



3M 



out 



(D.3) 



We may now show that Eq. (|D.3| ) is dynamically equivalent to the result (|2.26|) used in the 
text. We first integrate the second term in the integral above and obtain 



S 



(3) 
B 



8 7T 

— / " du \2R C 



3M out ] + — Rl Q ]n B 



M 



f 



R<? 



Mi 



(D.4) 



The last term does not affect the equations of motion, since the effective action is varied with 
fixed endpoints, and can therefore be dropped. It is shown in Section ^ that S^r depends on 
two arbitrary coefficients and for a = 1 and 7 = one just recovers the integral in Eq. ( p. 4 ). 



E Junction equations 

The junction equations at the shell surface S relate the jump in the extrinsic curvature Kij to 
the matter stress-energy tensor on S 1, 



[K 



out 
l 3Hn 



9ij S\ 



(E.l) 



where gij is the three-metric ( [2,5| ) in the proper time gauge N = 1 and the source term S. 



is given in Eq. Q2.7| ). Upon substituting the components (2.18) of the extrinsic curvature into 
Eq. (E.l) one obtains the 99- j unction equation 

16 7T 



E 



R[(3] 



in 
out ' 



(E.2) 



where we have introduced the shell proper mass E according to Eq. (|3.2| ). Analogously ( 2.19| ) 
yield the rr-j unction equation 



P 



2 

K 



(3 R M 
R + ^ + W]3 



+ 



M, 



out Pout 



00 



R 



out 



(3 out R-2M, 



(E.3) 



out 
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In Israel's coordinates (see Appendix |C|) the last term in the right hand side of Eq. (E.3) becomes 

UM out 1 V 2 dM out 



RP ou t RPout W dv 



(E.4) 

which renders our equation flE.3| ) the same as Eq. (3.3) of Ref. |13|] where the problem of inte- 
grating Eqs. ( |E.2| ) and ( |E.3| ) was addressed numerically for a particular choice of the luminosity. 

Let us focus on analytic results that can be obtained independently of the detail of the 
interaction between the shell and the out-flowing null dust. We take, for Qi n , flat Minkowski 
space (M in = 0) and 2M ou t(n) < Ri- 

1) Taking R > 2M out and R -C 1 in Eq. (|E.2j) gives, to leading order, E ~ (16 it/k) M outl which 
is the relation one would expect in (asymptotically) flat space. 

2) Expanding Eqs. (U) and (PH) for R - 2M out < 2M out gives, again to leading order, 



E~2M out (yi + R 2 + R) (E.5) 

— — 87rM out RP. E.6 

k R-2 M^ 

For E = it follows that M ou t = only if P always remains zero (dust). One can interpret this 
result as if, in general, the radiation can be fed by the gravitational energy of the shell whose 
extraction induces a surface tension (and changes the motion of R correspondingly). 

If we demand that P be regular everywhere, because of the denominator on the right hand 
side of Eq. QE.6|), for M out < the proper energy emitted per unit proper time diverges with 



R — > 2M ou t- This, together with Eq. (EJ3), namely E ~ M ou t, implies that E and M ou t would 
vanish at a value of R > 2M ou t. On the other hand, if the time derivative of the mass aspect 
vanishes according to the bound in Eq. ( 3.171 ) one (temporarily) recovers the non-radiating case, 



for which the shell crosses 2 M ou t at a finite proper time and with finite proper energy. The 
condition ( p. 17 ) is necessary for the stress-energy tensor of the radiation to be locally finite when 



the shell crosses the surface r = 2M out [31], that is 



\rprad\ I /\V I 

Hm , 1 uu 1 NO ~ lim , 1 ° utl - < oo . (E.7) 
i — >2M OU f (r — 2M ou t) Mout (R-2 Mout) 

We remark that the stronger condition M ou t = is necessary to ensure that r = 2 M ou t is a null 
surface after the shell has crossed it, which is a property of event horizons as opposite to apparent 



horizons [32]. In fact, the tangent to the surface r = 2M ou t in £7j n is = (Ti n , 2 M ou t, 0, 0) 
and has norm t^tu = — 1. In £l ou t the tangent would be f ~ (U, 2 M ou t, 0, 0) with norm 
tn ~ —4 M ou t U. Therefore f is time-like in Vti n and would become space-like (or null for 
2M out = 0) in n out ||. 

The mathematical origin of the bound ( |3. 17|) lies in the inadequacy of the coordinates (u, r) to 
cover both the interior and the exterior of r = 2 m in the Vaidya space-time (see Appendix |C|). It 
has also a physical interpretation in terms of the point of view of a static observer at r 3> 2 M OU f. 
On assuming that the metric for r S> 2M ou t can be approximated by the Schwarzschild line 

(E.8) 



element, the time too of the static observer is related to r according to [[y, 12 

dtoo (3 2 Mout R 



dr l-2Mout/R R-2M, 



out 
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Then Eq. (E.6) can be written as 



dE 32 ir dM out 

- — ~ — . (E.9) 

dtoo k dr 

This is in accordance with the fact that, due to the infinite redshift experienced by a distant ob- 
server, the luminosity of a star which collapses and forms a black hole would decay exponentially 
and eventually (for too +00) fade. On the other hand, if the flux of radiation measured by 
the distant observer does not vanish before all the proper energy of the shell has been radiated 
away, then the shell remains always outside the surface r = 2 M out until E vanishes (thus leading 
to flat Minkowski space instead of a black hole). 
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